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We study an optomechanical cavity, in which a buckled suspended beam serves as a mirror. The 
mechanical resonance frequency of the beam obtains a minimum value near the buckling temper- 
ature. Contrary to the common case, in which self-excited oscillations of the suspended mirror 
are optically induced by injecting blue detuned laser light, in our case self-excited oscillations are 
observed with red detuned light. These observations are attributed to a retarded thermal (i.e. bolo- 
metric) force acting on the buckled mirror in the inwards direction (i.e. towards to other mirror). 
With relatively high laser power other interesting effects are observed including period doubling of 
self-excited oscillations and intermode coupling. 

PACS numbers: 46.40.- f, 05.45.- a, 65.40.De, 62.40.+ i 



I. INTRODUCTION 

Systems combining mechanical elements in optical res- 
onance cavities are currently a subject of intense interest 
[IH^- Coupling of nanomechanical mirror resonators to 
optical modes of high-finesse cavities mediated by the 
radiation pressure has a promise of bringingthe mechan- 
ical resonators into the quantum realm @, |H 04HI (see 
Ref. fl4T | for a recent review). In addition, the micro- 
optoelectromechanical systems (MOEMS) are expected 
to play an increasing role in optical communications [l5j 
and other photonics applications [16l - ll8l |. 

Besides the radiation pressure, another important 
force that contributes to the optomechanical coupling in 
MOEMS is the bolometric force 0, H also 

known as the thermal force. This force can be attributed 
to the heat induced deformations of the micromechan- 
ical mirrors [2ll - |23l . [26j . In general, the thermal force 
plays an important role in relatively large mirrors, in 
which the thermal relaxation rate is comparable to the 
mechanical resonance frequency. Phenomena such as 
mode cooling and self-excited oscillations 
l3lj have been shown in systems in which this force is 
dominant 0, [H H H, |M, H . 

In this paper we investigate an optomechanical cav- 
ity having a suspended mirror in the shape of a trilayer 
metallic doubly clamped beam [34T - I361 ] . The system ex- 
perimentally exhibits some unusual behaviors. For exam- 
ple, contrary to other experiments, in which self-excited 
oscillations of the mechanical resonator are optically in- 
duced when the cavity is blue detuned [37| , here the same 
effect occurs for the case of red-detuned cavity. More- 
over, the dependence of the mechanical resonance fre- 
quency on laser power is found to be non-monotonic [38j . 
These findings are attributed to optically induced ther- 
mal strain in the beam, which gives rise to compressive 
stress that may result in buckling [39|-|44j. Bagheri et at. 
145] have recently reported the utilization of the buckling 
phenomenon to develop a non- volatile mechanical mem- 
ory element in a similar optomechanical system (see also 
Ref. H|). 

We generalized the theoretical model that has been 



developed in Ref. [32} to account for the effect of buck- 
ling [43, 13 • We find that close to buckling the effective 
thermal force acting on the beam can become very large, 
and consequently optomechanical coupling is greatly en- 
hanced. Consequently, the threshold laser power needed 
for optically driving self-excited oscillations in the mirror 
can be significantly reduced. We show that the proposed 
theoretical model can account for some of the experi- 
mental results. On the other hand, some other exper- 
imental observations remain elusive. For example, self- 
modulation of self-excited oscillations is observed in some 
ranges of operation. Further theoretical study, which will 
reveal the complex structure of stability zones and bifur- 
cations of the system [4j| , is needed in order to account 
for such findings. 



II. EXPERIMENTAL SETUP 

The experimental s etup (see Fig. [TJ is similar to the 
one employed in Ref. |33l. A fiber Bragg grating (FBG) 
[5fj| - l52| and a microlens |53j are attached to the end of 
a single mode fiber. The fiber can be accurately posi- 
tioned using piezomotors. The system, which consists of 
doubly clamped suspended multilayer micromechanical 
beam and the optical fiber monitoring assembly, is lo- 
cated inside a cryostat maintained under typical pressure 
of 10~ 3 mbar and temperature of 77 K. Optical cavity 
is formed between the freely suspended beam oscillating 
parallel to the optical axis of the cavity and between ei- 
ther the FBG or the glass- vacuum interface at the tip of 
the fiber, as shown in Fig. Ufa). 

The multilayer beam is fabricated by using bulk micro- 
machining process. A 200 nm thick Nb layer is coated on 
prefabricated silicon nitride membranes over silicon sub- 
strates using DC magnetron sputtering at working pres- 
sure of 5.2 x 10~ 3 torr and argon atmosphere. Pattern- 
ing on the coated substrates is done using photolithogra- 
phy followed by a liftoff process, in which a 10 nm thick 
layer of Cr and a 30 nm thick layer of gold-palladium 
(Auo.85Pdo.15) is patterned on top of the Nb layer. Front 
electron cyclotron resonance (ECR) plasma etching of the 
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FIG. 1: The experimental setup, (a) Sample and optical dis- 
placement sensor, (b) Side view sketch of the trilayer beam 
mirror, (c) Top view of the suspended trilayer beam A. 



unmasked Nb and silicon nitride is done to suspend the 
AuPd/Cr/Nb multilayer beam. Process is followed by a 
backside ECR plasma etching of the suspended beam to 
remove silicon nitride form back side. FigureQJb) depicts 
the sketch of the multilayer AuPd/Cr/Nb suspended mi- 
cromechanical beam having suspended length I = 222 /im 
, width b = 4/iin and total thickness h — 240 nm with 
effective Young's modulus E = 109.68 x 10 9 Pa, den- 
sity p = 9715.22 kg/ m 3 , thermal expansion coefficient 
a = 7.76 x lO^K -1 and Poisson's ratio a = 0.3961. 
These parameters are computed by averaging over the 
layers using expressions taken from Ref. 36]. The top 
view of the suspended trilayer beam mirror that is de- 
scribed in the following section is shown in Fig. [ljc). We 
will refer to this beam hereafter as beam A. 



The suspended multilayer beam can be capacitively 
actuated by applying voltage to it with respect to the 
ground plate on which the sample is mounted. The 
ground plate is located at 500 /im below the suspended 
beam. Incident optical power is controlled by attenuating 
the injected laser power of wavelength A ~ 1550 nm. The 
optical power reflected off the cavity is fed to a photode- 
tector, which is monitored using a spectrum analyzer, as 
shown in Fig. [Ha). 



III. SELF-EXCITED OSCILLATIONS INDUCED 
BY RED-DETUNED CAVITY 

As was pointed out above, the optomechanical cav- 
ity under study exhibits some unusual behaviors. As is 
demonstrated in Fig. [2] self-excited oscillations can be 
induced by red-detuned cavity. Subplot (a) shows the 
reflected optical power vs. the voltage V z , which is ap- 
plied to the piezomotor that is used to position the opti- 
cal fiber along the optical axis direction. The period of 
oscillation in the reflected power is A/2. In this experi- 
ment the wavelength A is not tuned to the reflective band 
formed by the FBG, and thus the optical cavity is formed 
between the freely suspended beam mirror and the glass- 
vacuum interface at the tip of the fiber. For this case the 
finesse of the cavity is much lower compared with values 
that can be achieved when the FBG is employed. Note 
that in addition to the oscillatory part, also the averaged 
(over one period of oscillation) measured reflected power 
exhibits dependence on cavity length [see subplot Ufa)]. 
The average value obtains a maximum when the cavity 
length coincides with the focal distance of the microlens. 
This occurs close to the value V z = 0. 

While the data seen in subplot HJa) is taken when the 
input optical power is attenuated well below the thresh- 
old of self-excited oscillations, in subplot d^b), which 
depicts the measured spectrum of reflected optical sig- 
nal, the input power is set to 0.16 mW. Self excited os- 
cillations of frequency of about 67 kHz are found near 
V z = 8 V and near V z = 28 V. In both these regions 
where self excited oscillations occur the reflectivity de- 
creases with V z , i.e. the cavity is effectively red-detuned. 
As was discussed above, the deviation from periodic de- 
pendence of the measured spectrum on cavity length with 
A/2 period is attributed to the effect of the microlens. 

The dependence on laser power and Pl is seen in Fig. 
|3l The threshold power of self-excited oscillations occurs 
at 0.092 mW. In subplot EJa) the static position yo of 
the center of beam A is measured vs. laser power below 
the threshold power. This is done by measuring the the 
reflected optical power vs. the voltage V z . For each value 
of the laser power this yield a plot similar to the one see 
in Fig. [2ja). However, the phase of oscillations is found 
to depend on laser power. The value of the static deflec- 
tion yo j which is extracted from that shift in the phase, 
is seen in subplot OHa) . The observed behavior indicates 
that the beam is deflected towards the optical fiber with 
increasing laser power. As will be seen below, this behav- 
ior indicates that the thermal force in the present case 
acts in the inwards direction, contrary to the more com- 
mon situation (e.g. for the case of radiation pressure), 
where cavity optomechanical forces act in the outwards 
direction, and lead to elongation of the cavity (rather 
than shortening as in the current case). 

In subplot [3jb) the mechanical resonance frequency /n 
is measured both below and above the threshold power 
of self-excited oscillations. While below the threshold 
of 0.092 mW the frequency /o is determined from forced 
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FIG. 2: (Color Online) Dependence on cavity length (beam 
A), (a) The photodetector DC voltage as a function of the 
voltage V z applied to the piezomotor holding the fiber. A 
change of 18.5 V in V z corresponds to a displacement of the 
fiber by A/2. The laser power in this measurement is set to 
the value of 0.08 mW, which is lower than the threshold of 
self-excited oscillations, (b) Spectrum analyzer measurement 
of the photodetector voltage measured with a laser power of 
0.16 mW. Self-excited oscillations are observed when the cav- 
ity is red-detuned. The white dotted line shows the calculated 
value of f2 c ff /2-7T for the regions in V z for which 7 e ff < 0. The 
following parameters have been used for the theoretical calcu- 
lation 9c = 0.06, 0t = 0.1, Q = 1.2 x 10" 4 , Cc = 2.8 x 10" 3 , 
7 = 0.5 x 10~ 9 , T B = 0.96, T A = 0.01, T R = 0.6, p Y = -1600 
and /5 T r = 2.6 x 10~ 7 . 



oscillations (FO) (i.e. from the peak in the measured 
frequency response), above threshold /o is determined 
from the peak in the spectrum of self-excited oscillations 
(SO). Note the non-monotonic dependence of fo on laser 
power. As will be argued below, this dependence suggests 
that the beam undergoes thermal buckling. 



IV. MECHANICAL EQUATIONS OF MOTION 

Most of the experimental findings that were presented 
in the previous section are not accountable by the theo- 
retical model that has been developed in Ref. |32[ (which, 
on the other hand, was very successful in accounting for 
previous measurements of Ref. [33j). We therefore gen- 
eralize the model to account for the effect of buckling in 
the mirror. While the derivation of the system's equa- 
tions of motion is described in the appendices, the final 
equations are presented below in this section. 

The height function y (x, t) is assumed to have the 
shape of the first buckling configuration of a doubly 
clamped beam [47], E8j. Thus, for the case where the 
clamping points are located at (x,y) = (±Z/2,0) the 
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FIG. 3: Dependence of laser power (beam A), (a) The static 
deflection j/o- (b) The mechanical frequency fo- The vertical 
red line represents the threshold of self-excited oscillations. 
Below the threshold fo is found from the peak amplitude of 
forced oscillations, whereas above the threshold fo is the fre- 
quency of self-excited oscillations. 



height function y (x,t) can be written as 

T 



l + cos • 



(1) 



where £ depends on time. The equation of motion for £, 
which is derived in appendix C, is given by 



£ + 2 7 £ - 
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where overdot denotes a derivative with respect to the 
dimensionless time r, which is related to the time t by 
the relation r = tj \fpA^/E, where p is the mass den- 
sity, A cs is the cross section area of the beam, E is 
the Young's modulus, 7 is the dimensionless damping 
constant, and F p and f2 p are the dimensionless ampli- 
tude and frequency respectively of an externally applied 
force. The equation of motion © contains two thermo 
optomechanical terms that depend on the temperature 
of the beam T. The first is the temperature depen- 
dent angular resonance frequency f2 and the second is 
the thermal force i*th- In terms of the dimensionless 
temperature 9 = (T — Tq) /To, where To is the temper- 
ature of the supporting substrate (i.e. the base temper- 
ature) the frequency f2 is given by f2 = f2o^(#), where 
fio = ypAcs/EoJo, wo is the mode's angular resonance 
frequency and the temperature dependence is expressed 
in terms of the dimensionless function v (9) . The di- 
mensionless thermal force is given by F t h = ^ 2 £c/y (9), 
where £c = f2o\/37u ( a ~ a s) £ 4 /167r 6 /, where a and a s 
are the thermal expansion coefficients of the metallic 
beam and of the substrate respectively, / is the moment 
of inertia corresponding to bending in the xy plane, and 
the dimensionless function fy (9) represents the beam's 
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temperature dependent deflection. For the case of a rect- 
angular cross section having width l y in the y direction 
and width l z in the z direction, and for the case of bend- 
ing in the xy plane I is given by I = lyl z /12. When 
small asymmetry is taken into account [see Eqs. (|C3[) 
and (|C5p ] the functions /y (9) and v (9) can be expressed 
as (see Fig. [5]) 



f Y (0) = Re (y/T^ T c 
f^(9) = -(9-9 c )+3.f Y (9) , 



(3) 
(4) 

where 9q is the dimensionless buckling temperature, 
which is related to the buckling temperature Tq by the 
relation 6c — (Pc — To) /To, and where the real dimen- 
sionless constant 9{ represents the effect of asymmetry. 



V. OPTICAL CAVITY 

The finesse of the optical cavity is limited by loss mech- 
anisms that give rise to optical energy leaking out of the 
cavity. The main escape routes are through the on fiber 
mirror (FBG or glass- vacuum interface), through absorp- 
tion by the metallic mirror, and through radiation; the 
corresponding transmission probabilities are respectively 
denoted by Tb, Ta and Pr. Let j/d be the displacement 
of the mirror relative to a point, at which the energy 
stored in the optical cavity in steady state obtains a lo- 
cal maximum. In the ideal case, all optical properties of 
the cavity are periodic in j/d with a period of A/2, where 
A is the optical wavelength (though, as can be seen in Fig. 
Ufa) , deviation from periodic behavior is experimentally 
observed) . 

It is assumed that j/d is related to £ by the following 
relation Anyn/X = (3y (£ + £f), where both f3y and £f 
are constants, which depend on the position of the fiber. 
For a fixed £ the cavity reflection probability Rc, i.e. 
the ratio between the reflected (outgoing) and injected 
(incoming) optical powers in the fiber, is given by (33j 

R ( T "- T r TR ) 2 + 2[l-cos(/3 Y (C + 6))] 

( TB+ r A+ T R ) 2 +2[1 _ CQs(M ^ + gf))] • 

The heating power Q due to optical absorption of the 
suspended micromechanical mirror can be expressed as 
Q = I (£) Pl , where Pl is the power of the monochro- 
matic laser light incident on the cavity and where the 
function I (£) is given by [33[ 

HC , , 2 -. (6) 



- Tb+Ta+T r 
2 



) + 2 [1- cos (/3 Y (£ + &))] 



VI. THERMAL BALANCE EQUATION 

The temperature T evolves according to the following 
thermal balance equation 

dP = Q _ H (T — Tp) 
dt ~ C C 



where C is heat capacity, H is thermal transfer coefficient 
and Q is the heating power due to optical absorption. In 
terms of the dimensionless temperature 9 and dimension- 
less time t the thermal balance equation becomes 



(8) 



where /3tr — y/T/EH/C is the dimensionless thermal 
rate, and j3p = •v/T / EPl/CT is the dimensionless in- 
jected laser power. 



VII. SMALL AMPLITUDE LIMIT 

The coupling between the equation of motion for £ [Eq. 
(|5J)] and the the one for 9 [Eq. (|SJ)] origins by three terms, 
the 9 dependent frequency f2, the 9 dependent force P t h 
(i.e. the thermal force) [see Eq. ([2])], and the £ dependent 
optical heating power [the term I (£) in Eq. {SJ] . The 
case where the first two coupling terms are linearized (i.e. 
approximated by a linear function of 9) is identical to the 
case that was studied in Ref. [32] , in which slow envelope 
evolution equations for the system were derived, and the 
amplitudes and the corresponding oscillation frequencies 
of different limit cycles were analyzed. By employing the 
same analysis (and the same simplifying assumptions) 
for the present case one finds that the coupling leads to 
renormalization of the mechanical damping rate, which 
effectively becomes 



7off = 7 + 



Q 2 



(9) 



where I' = d//d£ and where f Y = d/y/d^. The corre- 
sponding effective mechanical resonance frequency Sl c ff is 
given by 



a 



eff 



(7off - 7) 



TR 



fin 



(10) 



(7) 



The white dotted line in Fig. EJb) shows the calculated 
value of fieff [see Eq. ([TP)) ] for the regions in V z for which 
7cff < [see Eq. ©]. The system's parameters that 
were used for the calculation are listed in the caption of 
Fig. Note that Eqs. © and (JTDJ), which were de- 
rived by assuming the limit of small amplitudes, become 
inapplicable in most of the region, in which self-excited 
oscillations are observed in Fig. [2][b), and consequently, 
relatively large deviation between data and the predic- 
tion of Eq. (|10p is expected. 

The threshold of instability, i.e. Hopf bifurcation, oc- 
curs when 7 e ff vanishes. For the case where the dominant 
coupling mechanism between the mechanical resonator 
and the optical field arises due to radiation pressure, in- 
stability can occur only for the so-called case of blue de- 
tuning, i.e. the case where P < 0. However, contrary to 
the case of radiation pressure [37], HiJ , which always acts 
in the outwards direction, the thermal force due to buck- 
ling can act in both directions. We refer to the case where 
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FIG. 4: (Color Online) In addition to the harmonics in the 
measured spectrum due to self-excited oscillations of mode 1 
of beam B, subharmonics of order 1/2 and of order 1/5 are 
observed. 
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Laser Power [mVU] 

FIG. 5: (Color Online) The subharmonics of order 1/3 near 
input power of 11.9 mW are attributed to excitation of mode 
1 by the self-oscillating mode 2, which has a frequency fa that 
is 3 times larger than /i in that region. 



fy > as the case of outwards buckling, in which due to 
buckling the optical cavity becomes longer, whereas the 
opposite case of inwards buckling occurs when f Y < 0. 
As can be seen from Eq. ©, the threshold laser power 
/3p at which Hopf bifurcation occurs is inversely propor- 
tional to | /y | . Close to the buckling temperature 9q and 
when asymmetry is small (i.e. when |#f| <C 1) the factor 
|/ Y | can become very large [see Eq. ^ and Fig. [6], and 
consequently, the threshold laser power in that region can 
become very small. 



VIII. HIGH LASER POWER 

Other interesting effects have been observed with rel- 
atively high laser power. Two examples are presented 
below. The multilayer beam mirror that was used for 
the experiments that are introduced in this section (see 
Figs. |4] and [5J has the same layer structure as beam A 
[see Fig. [TJb)], however its length is / = 120 /an and its 
width is b = 20 /j,m. We will refer to this beam hereafter 
as beam B. The same fabrication process has been em- 
ployed for both beams. The two lowest lying modes of 
bean B, which are hereafter referred to as mode 1 and 
mode 2, have resonance frequencies fi and fa of about 
510 kHz and 1200 kHz respectively (note that these val- 
ues vary with laser power and cavity length). 

Self-excited oscillations of mode 1 are seen in Fig. [4] 
In this measurement the sample is kept grounded and 
the reflected optical power is measured using a spectrum 
analyzer as a function of input optical power. In the fre- 
quency window that has been employed for this measure- 
ment the first 5 harmonics of the self-oscillating mode 1 
are seen. In addition, for some values of input power the 
spectrum contains subharmonics. Close to input power 



of 11 mW the subharmonics indicate that period dou- 
bling occurs, whereas in the range 16.2 — 18.6 mW the 
period is multiplied by a factor of 5. While the period 
doubling effect has been predicted for a similar system 
by Blocher et al. [4^ | , the underlying mechanism respon- 
sible for the period multiplication by a factor of 5 is yet 
undetermined. 

The data in Fig. [5j which was taken with a modified 
cavity length, exhibits self-excited oscillations of mode 2 
in the laser input power range 6 — 16.5 mW. Above that 
range self-excited oscillations of mode 1 are observed. 
Near input power of 11.9 mW the ratio between and 
fi is 3, i.e. /2 = 3/i. This gives rise to strong intermode 
coupling, which leads to excitation of mode 1 by the self- 
oscillating mode 2, as can be seen from the subharmonics 
of order 1/3 that are measured in that range. 

Note that for other cavity lengths the measured spec- 
trum of reflected power can exhibit wide band response 
when the input laser power is sufficiently large. This pos- 
sibly indicates that the dynamics becomes chaotic [49| . 



IX. SUMMARY 

In summary, an optomechanical cavity having a buck- 
led mirror may exhibits some unusual behaviors when 
the operating temperature is close to the buckling tem- 
perature. In that range the instability threshold can be 
reached with relatively low laser power. Moreover, con- 
trary to the more common case, the instability is ob- 
tained with red detuned laser light for the case of inwards 
buckling. With relatively high laser power other interest- 
ing effects are observed in the unstable region including 
period doubling and intermode coupling. Further theo- 
retical work is needed to account for such effects. 
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Appendix A: Lagrangian 

Consider a beam made of a material having mass den- 
sity p and Young's modulus E. In the absent of tension 
the length of the beam is Iq. The beam is doubly clamped 
to a substrate at the points (x, y) = (±1/2, 0) and the mo- 
tion of the beam's axis, which is described by the height 
function y(x,t), is assumed to be exclusively in the xy 
plane. The corresponding moment of inertia is denoted 
by /. The Lagrangian functional C is given by [48[ 



1/2 

£= J dxL 

-1/2 

where 



(3 A El 



( 1/2 



V-v 



ox 



A csP fdy\ _N (dyY_EI (&V_ 
dt 2 \dx 2 \dx 2 



, (Al) 



-fy (A2) 



is the Lagrangian density, A cs is the cross section area of 
the beam and 



N = EA, 



l-lo 
lo 



(A3) 



is the tension in the beam for the straight case where 
y = 0. The beam's equation of motion is found from the 
principle of least action to be given by [13] 

^d 2 y l T A CS E W f , fdy\ 2 \ d 2 y 



-1/2 



d 4 n 



(A4) 



where T = pA cs is the mass density per unit length. 



Appendix B: First Buckling Configuration 



has 



Consider the case where the deflection y ( x, t 
the shape of the first buckling configuration |47|, i.e. 
y = yi (l + cos2y^), where y is a dimensionless time 
dependent amplitude. In order to account for a possible 
asymmetry, the force / is allowed to be a nonzero con- 
stant 148| ■ The Lagrangian (jAll) for the present case 
becomes 

3(3 A pl 5 (U) 2 



4tt 4 faEry* + fi z y - 



7T^ A Ei 3 y 4 



(Bl) 



where (3a = A cs /l 2 and where 0i 
£o can be expressed as 



I /I 4 . Alternatively, 



C a = T - U , (B2) 
where the kinetic energy To is given by 



Tn = 



(B3) 



where m = 3/?a/^ 3 /2 and where the potential energy Uo 
is given by 



where 



u = ^u(y), 



\y) = my + my 2 + ^y 4 



(B4) 



(B5) 



r]i = —2f/mouJQ, r)2 = 1 - /3l, m = ^ 4 /3AEl/m uj 2 , /3l = 
-N/4tt 2 ^E1 2 and where u 2 = (2tt) 4 ^E/3pA cs . 

Local minima points 34) of the potential Uq are found 
by solving = du/dy. The dimensionless potential u (y) 
can be expanded near one of its local minima points 34) 
to second order in y — 34) = £ as u = Uo + v 2 ^ 2 + O (£ 3 ) , 
where both uo and v = 7/2 + 6774^0 are constants. 

For the case where rjx = the coordinate 34) is given 
by 34) = ± R-c y/~ 772/2^4. The solution for the case of 
small 1 771 1 can be approximated by 



r\2 + i (^vlm) 
2 m 



1/3 



(B6) 



provided that 1 772 1 is small, i.e. provided that the sys- 
tem is close to buckling conditions. The correspond- 
in g approxim ation for v is found from the relation v — 
\J-q 2 + 6774^0 ■ 



Appendix C: Mechanical equation of motion 

Let a and a s be the thermal expansion coefficients 
of the metallic beam and of the substrate respectively. 
At some given temperature Ttf the tension N is as- 
sumed to vanish. Since the substrate is much larger 
than the beam, one may assume that the substrate ther- 
mally contracts (or expands) as if the suspended beam 
was not attached to it, namely the distance between both 
clamping points at temperature T = Ttf + At becomes 
l s = lo (1 + « s At), where lo is the distance at the ten- 
sion free temperature Ttf- The tension N needed to 
keep the beam clamped to the substrate is thus given 
by N — EA CS (a s - a) A T [55]. Thus, the dimension- 
less parameter /3l can be expressed as a function of the 
dimensionless temperature 



T-Tp 
To 



(CI) 
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where Tq is the temperature of the supporting sub- 
strate (i.e. the base temperature), as (3 L = 
PaPto (0 + Ptf) /4tt 2 /3i where /3 T o = T Q (a-a s ) and 
where (3tf = (To — Ttf) /To. Alternatively, /?l can be 
expressed as 



/3a/3to (g - 9 C ) 
4tt 2 /?i 



(C2) 



where #c = 47t 2 /?i//3a/3to — /?tf is the dimensionless tem- 
perature at which buckling occurs in the symmetric case 
(i.e. the case where 771 =0). 

With the help of Eq. (|C2p the parameters 34) an d 
v can be expressed as a function of the dimensionless 
temperature 0. For the case of small asymmetry one has 
[seeEq. jB6j)] 



= £c/y (0) 



(C3) 



where the function fy (9) is given by Eq. ([3]), 
4^/3 I (4 % 2 r, 4 ) 1/:; 
by 



' //3a/?to and the coefficient £c is given 



_ PaPto /3T (a - a s ) A cs pl 4 ul 

ec -y8^A^~V 16^7 • (C4) 

The corresponding approximation for v is given by 

v = vcfu (9) , (C5) 

where 



At: 2 Pi 



/3 a /3to _ /^cs/ 2 To (a - a.) 



4tt 2 / 



(C6) 



The 



and where the function f v (9) is given by Eq. (j4|. 
functions /y (0) and /„ (0) are plotted in Fig. [5] 

The temperature dependence of 3^o can be taken into 
account by adding a thermal force term Fth 0, HH, E3] 
to the equation of motion of the mechanical amplitude 
£ = y — yo- Similarly, the temperature dependence of 



v can be taken into account by including a thermal fre- 
quency shift term. When in addition damping and exter- 
nal driving are taken into account the equation of motion 
^ for £ is obtained. 
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FIG. 6: (Color Online) The functions f Y (6) and /„ (0) [see 
Eqs. (j3]) and ijljl] plotted for 3 different values of the asym- 
metry parameter Of and for the case where the dimensionless 
buckling temperature is 9c = 0.5. 
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